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It is exponentially hard to simulate quantum systems by classical algorithms, while quantum 
computer could in principle solve this problem polynomially. We demonstrate such an quantum- 
simulation algorithm on our NMR system to simulate an H2 molecule and calculate its ground-state 
energy. We utilize the NMR interferometry method to measure the phase shift and iterate the 
process to get a high precision. Finally we get 17 precise bits of the energy value, and we also 
analyze the source of the error in the simulation. 
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Introduction. - It is well known that quantum al- 
gorithms provides an excellent speedup towards classical 
algorithms in many problems [I] . Among these, the most 
famous example is the Shor's algorithm [2 , which solves 
the factoring problem exponentially faster than classi- 
cal methods. Besides the success of quantum compu- 
tation(QC) applied in computational problems, another 
important application of QC is the simulation of quan- 
tum systems, an idea conceived by Feynman[3 . Simula- 
tion of realistic quantum systems requires an exponential 
amount of resources to handle the large number of cou- 
pled Schodinger equations used to describe the particles 
in the system. However, within the QC architecture it 
was shown that universally simulating an physical sys- 
tem is efficient jl], i.e., costing polynomial resources with 
the size of the target system. 

Quantum simulation allows us to examine several 
new frameworks such as the many-particle system [5, 
pairing-Hamiltonian systems[7] and e£c.[5J [HI HO] . 
Specifically, several quantum-simulation algorithms were 
proposed to determine properties like thermal rate 
constant [H] and molecular energies [12] in quantum 
chemistry. The calculation of molecular energies is a 
fundamental problem focused in computational quantum 
chemistry. On classical computer, resources for a full sim- 
ulation of the molecular system scale exponentially with 
the number of atoms involved, limiting such full config- 
uration interaction calculations of molecular energies to 
diatomic and triatomic molecules [13 . However the cal- 
culation could be efficient if it is performed on a quantum 
simulator using the algorithm proposed by Aspuru-Guzik 
et al |12j . where information about the energy is trans- 
ferred onto the phase shift of a quantum register and 
measured by a quantum phase estimate algorithm (PEA) 
iteratively. 

Because of the practical importance of QS, experimen- 
tal demonstrations of these algorithms are not only of 
fundamental interests, but also a cornerstone of emerg- 
ing field of new technologies. Previous experiments 



have been performed to simulate some small systems on 
NMR|8, 9, 14 or ion-trap[15 platforms. However no such 
demonstrations have been performed till now due to the 
requirement of controlling sufficient number of qubits for 
simulation a molecule in original Aspuru-Guzik 's algo- 
rithm. 

In this letter, we adopt two methods to overcome the 
challenges: 1) choosing simplest case of the smallest 
molecule, i.e. H2] 2) utilizing the mature technology 
on NMR platform - the NMR interferometer - to mea- 
sure the phase shift. We also adopt a improved itera- 
tive scheme to extend the precision to a high level. By 
these methods, we realize this algorithm and calculate 
the ground-state energy of H 2 molecule to a precision 
of 17 bits. We also find that the precision of the result 
got from this algorithm is limited by the imperfection of 
simulating the operators in experiment. 

The Hamiltonian of H2 molecule. — Limited by the 
current technology, to calculate a large molecular energy 
in Full Configuration Interaction is not possible, so we 
choose the simplest situation: the ground-state energy 
of Hydrogen molecule in the minimal STO-3G basis [16]. 
The electron's Hamiltonian of H2 molecule with Born- 
Oppenheimcr approximation is shown as follows [TB] , 
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where Tj is the kinetic energy of the ith electron, and 
Vij is the coulomb potential energy between the ith elec- 
tron and the jth nucleus, while Oij is the coulomb po- 
tential energy between the zth and jth electron. In this 
two-nucleus and two-electron molecule in STO-3G ba- 
sis, each atom has a Is Gaussian-type function, and the 
two Is functions compose one bonding orbital with ger- 
ade symmetry and one antibonding orbital with ungerade 
symmetry. So there are 4 spin orbitals which can form 
6 configurations. Considering the singlet symmetry and 
the spatial symmetry of H 2 exact ground state, only two 
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configurations are acting in fact in the calculation: the 
ground state configuration |$o) an d the double excita- 
tion configuration l^j). Thus, the Hamiltonian matrix 
is(in atom units, the nucleus distance is 1.4a. u., only the 
electron's energy): 



H = 



22 1 
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0.1813 -0.2537,) 



(1) 



whose theoretical eigenvalue is -1.8516 a.M.. 

TTie calculation of molecular energy. — As shown in 
Figjlp, the calculation of the molecular energy in our ex- 
periment is described by four steps: 1) adiabatic prepa- 
ration of the system qubit to the ground state of the 
Hamiltonian H 2) application of the time evolution of 
the molecular Hamiltonian on the qubits to generate the 
phase shift on the probe qubit; 3) measurement of the 
phase shift on the probe qubit to extract the energy in- 
formation. We will introduce the three steps in detail as 
follows. 

In the first step, the quantum simulator's system qubit 
is prepared by Adiabatic State Preparation (ASP) to 
the ground state of the molecular Hamiltonian H. Ac- 
cording to the quantum adiabatic theorem \YI[ 118] . the 
qubit is prepared on an simple Hamiltonian's ground 
state and the systematic Hamiltonian of the qubit varies 
slowly enough from the simple Hamiltonian to H, if 
there's an energy gap between the ground state and the 
first excited state, the qubit will stay on the instanta- 
neous ground state of the system Hamiltonian. Thus the 
qubit finally is prepared on the ground state of H after 
the ASP. 

In the second step, An unitary operator U = e~ 
is applied to the state |^), with only generating a phase 
shift to the probe qubit by the controlled operation. Here 
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'' |*) where E = -2-k4>/t is the 



energy of H's ground state. Note that the energy E is 
negative so we make the phase </> to be positive and t is 
chosen properly to make the phase <f> ranges from to 1 . 

Finally in the measurement step, a four-bit inverse 
Quantum Fourier Transform (QFT) is adopted as the 
Relative Phase Measurement to evaluate the phase shift 
in the Aspuru-Guzik's proposal. This apparatus needs 
four qubits as probe qubits to get one precise bit with 
successful possibility of 15/16(1]. While in the NMR 
platform there's a mature technology NMR interferom- 
eter, named from the similar apparatus originally used 
in optics, which could easily measure the relative phase 
shift of the quantum states by modulating the spectrum 
patterns |19[ [20]. On our NMR interferometer, the phase 
shift could be evaluated with an error bound of less than 
±5°, much more precise than the performance of the orig- 
inal four-bit inverse QFT apparatus. Thus we utilize the 



interferometry to measure the phase shift in our experi- 
ment. 

The iterative scheme. — For useful practical applica- 
tion, it should be possible to iterate the above process 
to achieve arbitrary precision in the molecular energy. 
We made a small modification to the iterative scheme 
in Aspuru-Guzik's algorithm to improve its reliability 
As shown in Figjlji, for each iteration k we applied the 
controlled-J7fc and measure the phase shift. We start the 
iterations from Uq — U and iterate the process by choos- 
ing Uk+i = [e _ * 2ir ^ fc L/&] 2 ■ Here, n is the number of bit 
attained in each iteration and <p k = max{(j>k — 4> e rrbd, 0} 
where <f>k is the phase shift measured in the fcth itera- 
tion. Note that n is limited by the precision of the phase 
measurement in each iteration, i.e., 2~" > 2<^ erT .{ ) d- 
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FIG. 1: (a) General schematic circuit for calculating molec- 
ular energies, (b) Molecular structure of quantum register 
(CHCI3). (c) Pulse sequence for the adiabatic process to 
prepare the 13 C nucleus in the ground state of the molec- 
ular Hamiltonian H.(d) Pulse sequence to implement the 
controlled- Uk operation, where 9 = 0.226, 7 = 1.3458, /3 = 
0k ~ 4>k-i{4>k-i i 8 calculated from the result of last iteration 
and (3l = [2.0233,-2.6629,-2.4539,-0.7814,0.0322,0.2575]) 
for the kth iteration, and for the first control U operation, 
a — j and delay = jj — while a — 0,delay — for other 
5 times iteration. In fact, we could simplify the sequence to 
only three pulses on qubit 1 H except the first control U oper- 
ation. The U sequences are optimized as we chose a suitable 
r = 1.941122. 

The experimental implementation. - We used the 
13 C-labeled chloroform dissolved in aVacetone as a two- 
qubit NMR quantum computer, where 13 C nucleus was 
used for the system qubit while 1 H is for the probe qubit. 
The molecular structure is shown in FigJTJb). The natu- 
ral Hamiltonian of this two-qubit system is given by: 



it NMR = —cr. 
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where lo s /2tt and lu p /2tt are the Larmor frequencies of 
nucleus 13 C and 1 H, and J ps represents the J coupling 
constant, typically, J ps = 214.6iJz. The experiments 



3 



were carried out at room temperature on a Bruker AV- 
400 spectrometer. Let us now describe the experiment in 
detail. 

A) Preparation of the initial state: Starting from the 
thermal equilibrium state, we first created a pseudo-pure 
state (PPS) poo = + e| TTXTT I using the spatial 

average techniquc[21 , with I representing the 4x4 unity 
operator and e ~ 10~ 5 the polarization. However, the 
simulation algorithm requires the system qubit is in the 
ground state \ip g ) of the Hamiltonian of the H 2 molecule 
[Eq. (2)] while the probe bit is in the |+) = -±=| f) + | j). 

The |+) state was easily prepared by a pseudo-Hadamard 
gate Ry (ir/2) from | f) . On the other hand, the unknown 
ground state \ip g ) can be prepared by an adiabatic pro- 
cedure where we starts from a simple initial Hamilto- 
nian H whose ground state \ipo) is easily constructed. 
The adiabatic theorem tells that if the system Hamilto- 
nian H a d(t) varies slowly enough so that the adiabatic 
condition is satisfied |17j . the system remains in its in- 
stantaneous ground state. Therefore, if H a d(t) reaches 
the Hamiltonian H at t = T, the system is prepared 
the ground state \4> g ) of H. The time-dependent Hamil- 
tonian during the adiabatic passage is obtained by lin- 
ear interpolation: R a d = (1 — s)a x + sH with s = ^. 
Here we chose the initial Hamiltonian Ho — <J X whose 
ground state is |— ) = ^| f) — I I) prepared by the 
conjugated pseudo-Hadamard gate Ry(—n/2). We im- 
plemented the adiabatic preparation by discretizing the 
continuous adiabatic passage with the optimized param- 
eters [HI [53]: the discrete steps M + 1 = 6 and the 
total time T — 5.4ms. Therefore the unitary evolu- 
tion for the discrete adiabatic passage is then U a d = 

tM rrad T-fM „-iH m T 
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where the duration of each 



step is 5 = T/(M + 1). For each step, 
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7 *+0(6 3 ), (3) 

where s m = jj^T. Consequently, U%? can be im- 
plemented by the pulse sequence R_ x (9i) — R_ y (02) — 
R°{9 3 ) shown in Figjljc). The fidelity of the state ob- 
tained by the adiabatic preparation is around 98.93% 
comparing with the theoretical expectation 0.996. 

B) The controlled- [/fe operation. The controlled- Uk has 
the form of 



Z4 = | TXT \®U k . 



For the first iteration, 



k = 0, U = e 



-iHr 



Thus, the 



T 



Uo operation transforms the initial state tpi„ = 
) + | i))\t/> g ) into V/ = ^(| T) + e l2 ^| where 
4> = —Et/2tt with the energy E. As in an interferome- 
ter, the controlled logic gate effectively introduces a rel- 
ative phase shift 27r</> between 'two paths': the |0) and 
|1) states in the initially prepared superposition of the 
auxiliary qubit, which can be read out directly in NMR 
|19l I20j . The value of r, in principle, can be arbitrarily 



chosen to make <\> S (0, 1). For experimental convenience, 
we chose r = . = 1.941122, 

^/(2ff(l,2))2 + (fl-(l,l)-ir(2,2))2 

and the pulse sequence to implement the controlled- Uk 
operator is shown in Figjljd). The different Uk in each 
iteration is realized by adjusting the parameters a, (3 and 
delay . 

C) Measurement. The relative phase shift is obtained if 
we measure the NMR signal of the auxiliary qubit ((c~)): 

far) = (V'/ko 1^/) = cos(2tt0) +isin(27T(/>). 

As a result, the quadrature detection in NMR serves as 
a phase detector, i.e., the Fourier-transformed spectrum 
gives the relative phase information. Here we take the 
initial state of ipi n as the reference phase. The experi- 
mental spectra for each interaction are showed in Fig(2] 
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FIG. 2: Experimental 1 H spectra for all six iterations along 
with the reference spectrum of the initial state (denoted 
by k — —1 as the dashed line). The following spectra 
(k = ~ 5) are observed after each iteration, which pro- 
vides the relative phase information. The phase measured in 
this iteration is used in the next iteration. 

After each time iteration of the above procedure, we 
measure the phase shift and prepare the operator for the 
next iteration. After measured the phases, we uses a 
recursive method to rebuild the <f> as the experiment re- 
sult. The recursive method is formulated as follows ,for 
the result 4> exp after /cth iteration, 



rbd - 



4>1 is the intermediate value only for calculation. The 
recursive variable i iterates from i — k to 1 with = <pk 
where <j>k is the measured phase in the /cth iteration. And 
finally we get the result by <j) e xp — 4>o- 

The result of the iteration is shown in Fig(3^i and Tab|T] 
The value ifiexp w e got from experiment approaches the 
theoretical value of 4>th rapidly as the iterations proceeds. 



However, there's still a gap between 
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FIG. 3: (a) The energy values got from the experiment in 
each iteration. The green line is the spline fit of experimental 
values and the dash red line is theoretical expectation. The 
experimental value approaches theoretical expectation expo- 
nentially. The inside figure shows the small gap between ex- 
perimental value and theoretical expectation, (b) The error of 
phase shift generated by Uk in the iterations. The error scales 
exponentially with k due to the imperfection of the operator 
U. 



TABLE I: The Theoretical Expectation and Experimental result 
of <j> after each Iteration. The experiments get 3 bits of <j> in each 
iteration. The bits with an underline in each <p ex p denotes the 
precision of the value. Overall, we get 17 bits of precision by 6 
iterations. 
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is larger than 2 -18 after 6 iterations. This means that 
although we get 18 bits from the experiment, only 17 bits 
of them are precise. 

This is because the time evolution of molecular Hamil- 
tonian H is simulated by U — e~ lHr in the experiment 
and this simulation could not be arbitrarily precise. Thus 
the error of Uf., which contains 2 k 's power of U, scales 
exponentially with k. We figure out this error in Fig|3]3, 
which shows that in the 5th iteration the error is large 
enough and we can not get any more precise bits. Fi- 
nally the molecular energy got from our experiment is 
-1.851569, with 17 precise bits towards the theoretical 
value -1.8516. 

Conclusions. — We have demonstrated a quantum al- 
gorithm to calculate the molecular energies on an 2-bit 
quantum computer. This is one of early researches of 
quantum-simulation algorithms for chemical-interested 



problems and could move this field significantly forward. 
We also found that the precision of interested properties 
attained by this algorithm relies on the precision of imple- 
menting the time evolution of the molecular Hamiltonian 
in the experiment. Anyway, we made this early exper- 
imental progress towards the long-term goal of a rapid 
quantum chemistry calculation by quantum computers. 

Note that, we have learned of a similar work [24 done 
concurrently on photonic system since writing this paper. 
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